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ABSTRACT 

Context. Since the discovery of its dusty disk in 1984, /3 Pictoris has become the prototype of young early-type planetary systems, and 
there are now various indications that a massive Jovian planet is orbiting the star at ~ 10 au. However, no planets have been detected 
around this star so far. 

Aims. Our goal was to investigate the close environement of /3 Pic, searching for planetary companion(s). 
Methods. Deep adaptive-optics L'-band images of p Pic were recorded using the NaCo instrument at the Very Large Telescope. 
Results. A faint point-like signal is detected at a projected distance of ^ 8 au from the star, within the North-East side of the dust disk. 
Various tests were made to rule out with a good confidence level possible instrumental or atmospheric artefacts. The probability of a 
foreground or background contaminant is extremely low, based in addition on the analysis of previous deep HST images. Its U = 1 1.2 
apparent magnitude would indicate a typical temperature of ~ 1500 K and a mass of ~ 8 Mjup. If confirmed, it could explain the main 
morphological and dynamical peculiarities of the Pic system. The present detection is unique among A- stars by the proximity of the 
resolved planet to its parent star. Its closeness and location inside the j3 Pic disk suggest a formation process by core accretion or disk 
instabilities rather than binary like formation processes. 
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1 1. Introduction 

' Understanding planetary systems formation and evolution has 
become one of the biggest challenges of astronomy, since 
I the imaging of a debris disk around JS Pictoris in the 80's 
. (ISmith & Terrile 19E4^ and the discovery of the first exo- 
" planet around the sola r-like star 51 Pegasi during the 90's 
( Mayor & Queloz 1995| ). While about 20 debris disks - disks 
containing dust which is not primordial but produced by colli- 
sions among larger rocky bodies - have been resolved at optical 
wavelengths today, Pic, a A5V star at a distance of 19.3+0.2 pc 
(ICrifo et al. 1997K remains the best studied young (12^^ Myr; 
IZuckerman et al. 2001 ) system, with an impressive amount of 
indirect signs pointing toward the presence of planets. 

The disk shows a relative inner void of matter inside 
50 astronomical units (au). Lecavelier des Etangs et al. 
(1995) presented intriguing light variations possibly due to 
disk inhomogeneities produced by a Jupiter size planet at 
> 6 AU. Several asymmetries have been identified in the 
disk at optical (IKalas & Jewitt 19951 [Heap et al. 2000| ) and in- 

Send offprint requests to: A.-M. Lagrange 

* Based on observations collected at the European Southern 
Observatory, Chile, ESO; runs 072.C-0624(B) and 60.A-9026(A); and 
observations made with the Space Telescope Imaging Spectrograph on 
board the NASA/ESA Hubble Space Telescope. 



frared (Telesco et al. 2005) wavelengths, as well as a warp at 
~ 50 AU (Mouilletetal. 1997; Heap et al. 2000| . The struc- 
ture is well reproduced by the deformation induced on col- 
liding planetesimals by a giant planet on a slightly in- 
clined orbit within 50 au from the star (IKristetal. 1996[ 
pouillet et al. 1997 Gorka vyi et al. 2000[ [Augereau et al. 2001 
and T hebault & Beust 200 iT Silicate dust is observed as 
circumstellar rings at 6, 16, and 30 au from the star 
(lOkamoto et al. 20041) , which could be explained by the pres- 
ence of a 2-5 -Jovian-mass (Mjup) planet at ~ 10 au from the star 
( Freistetter et al. 20071) . Evaporating star grazing comets have 
been evidenced (see, e.g., [Lagrange et al. 2000| for a review) and 
dynamical simulations showed that the gravitational perturba- 
tion of at least one giant planet at ~ 10 au can account for the 
observed rate of evaporating bodies (^ Beust & Morbidel li 2000b . 
However, no planets have been detected so far, either through 
direct imaging or through radial velocity studies, due to the in- 
strumental limitations of both techniques dGalland et al. 2006b . 
In particular, the high spatial resolution imaging detection capa- 
bilities were so far limited to distances typically ^ 15-20 au. 
We used the NAOS-CONICA instrument (NaCo), installed 
on the Very Large Telescope UT4 (Yepun) set in Paranal 
(Chile), to benefit from both the high image quality pro- 
vided by the Nasmyth Adaptive Optics System (NAOS; 
.Rousset et al. 2003i) at infrared wavelengths and the good dy- 
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namics offered by the Near-Infrared Imager and Spectrograph 
(CONIC A; Lenzen et al. 2QQ3I) detector, in order to study the im- 
mediate circumstellar environement of p Pic. 



2. Observations and data reduction procedures 

2.1. Observations 

L'-band images of p Pic (y = 3.8, L'=3.5) were obtained be- 
tween 2003 November 10 and 2003 November 17 with NaCo. 
The visible wavefront sensor was used with the 14 x 14 lenslet 
array, together with the visible dichroic. We used the CONIC A 
L27 camera, which provides a pixel scale of ~ 27 mas. Saturated 
images of p Pic were recorded, with detector integration times 
(DITs) of 0.175 s and number of detector integrations (NDIT) of 
100 or 200. Every two exposure^, spatial offsets were applied in 
order to allow sky and instrumental background removal. Non- 
saturated images were also recorded to get images of the stellar 
point spread function (PSF) as well as a photometric calibration. 
In such case, we added the Long Neutral Density filter (transmis- 
sion ~ 0.018) in the CONIC A optical path, and recorded images 
with DITs of 0.4s. 

In addition, the binary IDS 22141S3712 (separation 
p = 6 630 + 10 mas, position ang le PA = 302.06 ± 
0.07°; van Dessel & Sinachopoulos 1993|| ) was observed on 
November 1 1 as an astrometric calibrator. A mean plate scale of 
27. 105 + 0.041 mas and a true North orientation of -0. 10 ± 0.07° 
were derived and used to calibrate all p Pic images. Saturated 
and non saturated exposures were taken on the reference star 
HR 2435 (AOII, V = 4.4). The purpose is to correct for the star 
halo (the wings of the PSF) present in the saturated exposures. 
To optimize the removal of any fixed speckle, Pic and HR 2435 
were observed, as much as possible, at close parallactic angles. 
Finally, twilight flat fields were also recorded in U band. 

The observing conditions varied from exceptional (coherent 
energ}Q EC > 70%, coherent time tq > 20 ms) to reasonable 
(EC ~ 50%, To of a few ms), and sometimes poor (EC ~ 20- 
35%, To ~ 1-2 ms), over the run; the data quality varies accord- 
ingly. The best data set on p Pic was obtained on November 10. 
In the following, we will describe the data reduction and analysis 
of three sets of data: (A) the very best set obtained on November 
10; (B) a set of data with slightly poorer image quality obtained 
on the same night, with a shorter total exposure time, and (C) a 
set of data with poorer image quality obtained on November 13. 
However, Set C is representative of the best data obtained in the 
nights following November 10. The instrumental configurations 
and the status of image quality for these three data sets can be 
found in Table [H 

2.2. Data processing 

The first step was the cosmetic correction (bad-pixels, flat- 
fielding, background subtraction) and recentering of individ- 
ual offset positions of p Pic and HR2435 observations. A 
first method was to directly apply standard routines from 
the eclipse library (IDevillard 19971 ), using classical cross- 
correlation algorithm. A special care was taken to estimate the 
background at each given position by averaging images obtained 




^ An exposure is completed after DIT x NDIT. 

^ Since it is not possible to measure the Strehl ratio on our saturated 
data, we can only use the information provided by NAOS to assess the 
image quality. The coherent energies are those measured by the system 
in K band. 



Fig. 1. p Pic and HR2435 recentered and saturated U images 
(top left and top right, respectively) in data set A. Below are 
the divided (bottom left) and subtracted (bottom right) images. 
North is up and East is to the left. A candidate companion is 
clearly detected at a PA of ^ 32°, i.e., along the NE side of the 
disk, at a separation of about 0^.^41 from the star. 



at the previous and successive offset position in the observing 
sequence. Alternatively, a second method was used with an im- 
proved software dedicated to adaptive optics (AO) image pro- 
cessing (see Gratadour et al. 20 05 ). Following a different ap- 
proach for the background subtraction, images at individual off- 
set positions were recentered using a maximum likelihood al- 
gorithm at the level of the tenth of a pixel or better. The same 
overall process was used for both the object and the reference. 
Consistent results are found in terms of recentering precision 
and background subtraction (at less than the backgroung noise 
ofO.9 ADU). 

As a second step, three parallel approaches were followed to 
study the close environement of p Pic: 

- first approach consists in removing the PSF wings from the 
saturated images of p Pic by a simple minimization of the 
residuals. To do so, we first divided p Pic images by the 
ones of HR 2435 obtained under similar conditions (same 
DIT, NDIT, and pupil position). We then computed the scal- 
ing factor to be applied to the reference in order to scale its 
flux to that of p Pic. We then subtracted the reference images 
to the p Pic ones. Recentering and scaling processes are re- 
peated to minimize the residuals with respective precisions 
of sub-pixels and 5%. Tests have been performed using ro- 
tated images of p Pic for PSF subtraction. However, the fact 
that PSF halo is not centro- symmetric and that the statical 
aberrations are not overlaid worsen the subtraction process. 

- the second approach is to follow the same subtraction se- 
quence but applying the maximum likelihood algorithm of 
IGratadour et al. 2005l (see Fig.[T]). The algorithm is used here 
to recenter the reference star with p Pic and therefore confirm 
the previous estimation. Similar precisions are achieved. 

- the last approach was actually to use the MISTRAL decon- 
volution algorithm ( [Mugnier et al. 2004] ), based on a maxi- 
mum a posteriori scheme. Nevertheless, MISTRAL relies on 
a strict convolution process between image and reference 
which is not the case for our saturated data. A first step is 
therefore to perform a posteriori correction of saturated part 
of the image and reference. This is done using a simulated 
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COMPANION 



Set 


Star 


Date 


DIT 


NDIT 


^exp 




secz 


{EC)^ 


(to)I 


Separation 


PA 




AL' 








(s) 




(s) 


(0) 




(%) " 


(ms) 


(mas) 


(0) 




(mag) 


A 


ygPic 


2003-11-10 


0.175 


200 


665 


-18/-9 


1.125/1.119 


62.71 


11.47 


411 + 8 


31.8 


+ 1.3 


7.7 + 0.3 




HR 2435 


2003-11-10 


0.175 


200 


630 


-20/- 12 


1.148/1.140 


63.06 


10.65 










B 


/?Pic 


2003-11-10 


0.175 


100 


420 


-27/- 19 


1.138/1.127 


57.94 


8.08 


411 + 8 


31.5 


+ 1.3 


7.9 + 0.4 




HR 2435 


2003-11-10 


0.175 


100 


420 


-28/-20 


1.161/1.149 


60.23 


9.70 










C 


pVic 


2003-11-13 


0.175 


200 


385 


-24/- 19 


1.161/1.149 


46.94 


4.50 


401 + 8 


32.1 


+ 1.4 


7.6 + 0.4 




HR 2435 


2003-11-13 


0.175 


200 


350 


-22/- 17 


1.150/1.144 


44.99 


2.72 











^ Range of parallactic angles at the start/end of the observation. 
^ The average coherent energy as estimated on-line by AO. 
^ The average coherence time as estimated on-line by AO. 




Airy pattern. The top of the Airy pattern replaces the im- 
age saturated pixels. The flux level is adjusted using the first 
Airy rings. Such a correction is possible because of the very 
good Strehl ratio on the image. Meanwhile, if this a poste- 
riori correction does not significantly aff'ect the restitution 
of the object structures, it could obviously degrade the rela- 
tive photometry. Hence, the deconvolution process is an al- 
ternative approach (compared to reference subtraction and 
division) for the image processing which is less sensitive to 
image versus reference centering. It allows us to provide the 
best measurements of relative position with a precision of 
0.3 pixels, but remain uncertain for the relative photometry 
due to the use of saturated images. 




3. A candidate companion in the jS Pic disl<? 

3.1. Results from the highest-quality data (set A) 

Using three independent approaches, the companion detection is 
confirmed at the same location. The resulting images, using the 
maximum likelihood algorithm for recentering, are reported in 
Figure[T] The companion candidate (hereafter, the CC) point-like 
signal is clearly visible in the divided and subtracted images. The 
maximum of the signal is about 190 ADU. Diff'erent techniques 
(variable aperture photometry, 2D-gaussian fitting, PSF-fitting) 
were used to extract the CC flux, giving consistent results. As the 
reference recentering and rescaling actually dominate the flux 
measurement precision, flux uncertainties were derived consid- 
ering respective variations of 0.3 pixel and 5% in the subtraction 
process. We obtain a contrast of AU = 7.7 ± 0.3 between the 
CC and j3 Pic. Using deconvolution, we derive a separation of 
411 + 8 mas and a PA of 31. 8 + 1.3° relative to the primary, i.e 
along the NE side of the disk. 

The use of difl'erent methods excludes artefacts created dur- 
ing the reduction process. In particular, the result of the deconvo- 
lution rules out any efl'ect that could be introduced by imperfect 
estimation of the ofl'set between jS Pic and HR 2435 saturated 
images due, for instance, to a possible contribution of the disk. 
We did check anyway that the disk signal is very faint and not 
significantly asymmetric. To rule out detector eff'ects, we looked 
for possible remanence and electronic ghosts that could occur 
because our images are saturated. Individual images inspection 
exclude any contamination by these two efl'ects that rapidly dis- 
appear after a few frames. Artefacts due to the very good but 
still imperfect AO correction are still possible. However, aber- 
rations due to a modulation of the deformable mirror would 
generally lead to the presence of either symmetrical or anti- 
symmetrical patterns. Static aberrations should be present equiv- 




Fig.2. Top left: simulated planets at PA of 150°, 210° and 
330°. Top right: composite image of ft Pic plus the fake planets. 
Bottom left: division of the composite image by the saturated 
image of HR 2435. Bottom right: scaled subtraction of the com- 
posite image by the saturated image of HR 2435. Note that even 
a slight (0.3 pixel) relative ofl'set between JS Pic and HR 2435 
impacts the resulting shape of the fake planets as much as the 
candidate one. In particular, triangular shapes can be observed, 
due to the proximity of the slightly inner Airy ring. 



alently around /3 Pic and HR 2435, as both stars were observed 
with similar pupil configurations. We have then tested the possi- 
ble impact of an imperfect removal of static speckles due to the 
variation of the parallactic angle during the observations of j3 Pic 
and HR 2435 (up to 8°). We processed individual pairs of data 
of JS Pic and HR 2435 taken with parallactic angles equal within 
+0.4° and added up the individual subtracted images. The CC 
is still present and appears slightly sharper (but still compatible 
with the instrumental resolution). In addition, since the same sig- 
nal is also observed a few nights apart (see below), we conclude 
that quasi- static aberrations are unlikely. 

To further assess the reality of the detection and test the CC 
photometry, we added three 'fake planets' at similar separations 
but diff'erent PA (150°, 210°, and 350°) to the recentered and 
stacked image of JS Pic. The fake planet images were generated 
by scaling and shifting an unsaturated image of Pic taken dur- 
ing the same night. To match the level of the observed signal, the 
magnitudes of the fake planets are scaled to the measured flux 
ratio on set A. We then subtracted the scaled image of HR 2435 
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Fig. 3. Top: Scaled subtraction of the composite image by the 
saturated image of HR 2435 for sets A, B and C. Bottom: 
Simulated planets at PA of 150°, 210° and 330°.Note that even 
a slight (0.3 pixel) relative offset between p Pic and HR 2435 
impacts the resulting shape of the fake planets as much as the 
candidate one. In particular, triangular shapes can be observed, 
due to the proximity of the slightly inner Airy ring. 

to that of the composite image. The result is shown in Figure O 
from which it is clear that the fake planets produce features sim- 
ilar to the observed signal, supporting our contrast estimate. We 
therefore, derive an apparente magnitude of = 11.2 + 0.3 for 
the CC. 

3.2. Results from lower quality data (sets B and C) 

All other sets of data were processed as before. It quickly ap- 
peared that the detection becomes more marginal when the expo- 
sure time and/or the image quality decrease. For the best remain- 
ing sets in terms of atmospheric conditions and relative parallac- 
tic angle between JS Pic and HR 2435 (described in Section |2J1 
and Table [TJ, the results are shown in Figure [3] The same point- 
like signal is seen in both cases, with however a lower signal- 
to-noise ratio. This is due to the lower exposure times and the 
(slightly) poorer image qualities. This is illustrated in Figure O 
where we also added fake planets with a magnitudes equal to that 
of the fake planets used for Set A to the data sets B and C. The 
fake planet signals are, as the CC signal, much less detectable 
under these conditions; hence the results from our 3 data sets are 
consistent. 



4. Bound companion or background object? 

With these data alone, it is not excluded that the CC could be a 
foreground or background object. To definitely test this hypoth- 
esis, second epoch measurements are needed. Using the known 
13 Pic proper and parallactic motions, we have considered the po- 
sition of the CC assuming it would be a stationary contaminant. 
It would be closer than 0'.'2 to Pic in 2008 and hence, not de- 
tectable with current technics. In past years, j3 Pic has been mon- 
itored by numerous programs of the Hubble Space Telescope, 
including observations using the Advanced Camera for Surveys 
(Golimowski et al. 2006 and Kalas et al. 2005), the Near-Infrared 
Camera and Multi-Object Spectrometer (Brown et al. 1999) and 
the Wide Field Planetary Cameras (Kalas et al. 2000). The lack 



of spatial resolution or the size of the coronographic mask pre- 
vented from getting any reliable hints of any point sources close 
to the star. 

To our knowledge, the highest angular resolution and dy- 
namical data available are the coronagraphic data taken in 1997 
September with STIS (|Heap et al. 2000). These authors were 
able to probe the close stellar environment down to 0'.'75 from 
JS Pic. They did not report any CC in their images at the expected 
location of a.'9 North and a.'25 East to JS Pic. Based on STIS de- 
tection limits and the brightness derived by Heap et al. (2000), 
all objects brighter V < 17 would have been detected. Red gi- 
ants or supergiants can be ruled out for their distances that would 
bring them unrealistically far (typically ^ 10 000 pc) from us. If 
we consider a foreground or background stellar field contami- 
nant, it would therefore have a large (V - L > 5.8) color, i.e a 
spectral type later than mid-M. Based on a number density of L 
dwarfs of 1.9 x 10"^ pc"^ given by Burgasser (2001) and Cruz et 
al. (2003), the probability of finding a foreground or background 
field L dwarf in a region of 500 mas radius around jS Pic and with 
L' = 11.2 is about 10"^^. Without any assumption on the con- 
taminant spectral type, one can use galactic population model 
outputs to estimate the probability to find any (V < 12) galac- 
tic source at 500 mas from JS Pic. We find a low probability of 
6 X 10"^. Last, we cannot strictly rule out a contamination by an 
extragalactic source, such as a high-redshift quasar. In conclu- 
sion, a contamination appears very unlikely. In addition, the fact 
that the candidate companion falls into the disk strongly favors 
that it is bound to the star. 



5. Implications on the understanding of the ^ Pic 
system 

With an U magnitude of 1 1.2 ± 0.3, and assuming a distance of 
19.3 ± 0.2 pc and an age of 12+^ Myr (Zuckerma n et al. 20011) . 
the mass of the CC is estimated to 9^2 ^Jup using COND mod- 
els (Baraff'e e taL200 3| and 8^^ Mj^p using DUSTY models 
(Chabrier et al . 20001) . The planet should still be in the phase of 
cooling: Dusty and COND models predict eff'ective temperatures 
of 1 400 and 1 600 K, respectively. The validity of these models 
in the case of young planets formed by core accretion has been 
recently questionned (Marl ey et al. 2007) ) on the basis that the 
accretion shock might impact the planet initial internal entropy 
and its subsequent early thermal evolution. These authors claim 
that young giant planets could be significantly cooler and fainter 
than predicted so far. However, the treatment of the accretion 
shock is still a matter of debate. The impact of the initial disk 
mass is also to be studied. This may be important in the present 
case as j3 Pic is significantly more massive than the Sun. In the 
present case, |Marley et al. 2007| model predicts a luminosity ten 
times fainter from what we observed for an 8 Mjup planet at the 
age of /3 Pic. However, we can already note that a companion 
at a true separation of ~ 8 au could not be significantly more 
massive than 10-20 Mjup since otherwise it would have been 
detected through radial velocity measurements, as shown by a 
detailed analysis that will be presented in a forthcoming paper 
(Desort et al., in preparation). 

If the observed projected separation happens to be the phys- 
ical one or close to the physical one, then the companion ex- 
plains three very intriguing and so far unique characteristics of 
the p Pic system: the warp in the inner disk, the inner belts and 
the falling evaporating bodies (IFreistetter et al. 20071) . Mouillet 
et al. (1997) and Heap et al. (2000) showed that the warp of the 



A.-M. Lagrange et aL: A probable giant planet imaged in the p Pictoris disk 



5 



disk constrains the (mass M^, semi-major axis a) domain of the 
planet as follows: 

log(M^/M^) + 21oga + logr ^ 6.7, (1) 

where = 1.8 M© dCrifo et al. 1 997) and t is the age of the 
star. Given an age ^ = 12 Myr, masses Mp = 6 and 13 Mjup give 
a = 9.7 and 7.6 au, respectively, nicely bracketting the measured 
projected separation. This further strengthens the likely hood that 
the observed CC is the planet causing the warp. 

At the time of the submission of this letter, this companion 
was to our knowledge the first one possibly detected around an A 
type MS star. Recent results published since then by Marois et al. 
(2008) and Kalas et al. (2008) show also the presence of planets 
at distances 24 to 118 AUs from two A types stars. This com- 
panion could be anyway the first extrasolar planet ever imaged 
so close to its parent star: 8 au. In particular, it would be located 
well inside the orbits of the outer planets of the Solar System. 
Its closeness and location inside the JS Pic disk suggest a for- 
mation process by core accretion or disk instabilities rather than 
binary like formation mechanisms proposed for the companions 
to 2M1207 and AB Pic (Chauvin et al. 2005a, b). Further direct 
imaging observations should allow to constrain the planet orbit 
and hence measure the planet dynamical mass. This would in 
turn allow to test the models of planet formation. This is dramat- 
ically needed as the evolutionary models are not calibrated with 
real data in these ranges of masses and ages. Finally, we note that 
unfortunately, if the projected separation is the physical one, the 
planet orbital period is about 16 years and should be at present 
time much closer to the star and hence undetectable with NaCo. 
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